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Abstract 

We calculate the rotation of the inertial frame within an almost flat cylindrical region 
surrounded by a pulse of non-axially-symmetric gravitational waves that rotate about 



^vQ the axis of our cylindrical polar coordinates. Our spacctimc has only one Killing vector. 

> 

r- — It is along the z-axis and hypcrsurfacc orthogonal. We solve the Einstein equations to 

first order in the wave amplitude and superpose such linearized solutions to form a wave 
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pulse. We then solve the relevant Einstein equation to second order in the amplitude to 
find the rotation of inertial frames produced by the pulse. The rotation is without time 



The influence of gravitational wave angular momentum on the inertial frame demon- 

X 

!-h strates that Mach's principle can not be expressed in terms of the influence of the stress- 



energy- momentum tensor alone but must involve also influences of gravitational wave en- 
ergy and angular momentum. 
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1 Introduction 



Mach [21] rejected Newton's [26 absolute space and suggested that some average of the 
positions and motions of the masses in the Universe determines the inertial frame. 

Einstein [ID], jllj was inspired by Mach's principle and his earlier work [9J showed that 
the local inertial frame was not absolute but was affected by the rotation and acceleration 
of matter nearby. Extrapolating from this he said, in essence, the relativity of inertia is not 
fulfilled if inertia is affected by the presence of matter it must be entirely caused by it. He 
also said there is no inertia of mass against space but only inertia of mass against mass. Thus 
he hoped that the inertial frame might be determined by the distribution of the stress tensoiQ 
Tuv- His efforts to remove the external influence of the boundary conditions at infinity on the 
solutions of his equations led him to consider closed spaces with no spatial boundaries |10] . 
Here we show that within General Relativity, any general statement of Mach's principle that 
attributes all dragging of inertial frames solely to the distribution of T^ v as the origin of inertia 
is false (see [5]). One counter example suffices to demonstrate this. We find that there is 
an almost flat cylindrical region near the axis of a revolving gravitational wave pulse (which 
inevitably has no T^ v ) and demonstrate that the inertial frame in the cylindrical interior 
rotates relative to the inertial frame at great distances. 

The exact cylindrical wave solutions [2] , [12] , [7J , [30] carry no angular momentum about 
the axis, because their wave fronts have normals with no component around the axis. Bondi 
[6] emphasises this point. 

In a sense our work pushes at an open door, since now all accept that gravitational 
waves exist and carry both energy and angular momentum. This is not merely known theo- 
retically, but also observationally [28], [27], cf. also [M] and [20 . 

It would be astonishing if that energy and angular momentum did not produce gravity 
and rotation of the inertial frame. Indeed, as we described in more detail in paper [1] , the 
work of Corvino [8] demonstrates that this must be the case and in [1] we have explicitly 
calculated the gravity of a non- linear pulse of gravitational waves. 

Here our aim is to produce a nice clean example of the rotation of the inertial frame in 
an almost flat region surrounded by rotating gravitational waves. 



1 Indices X, jj,,u, p, ■ ■ ■ = 0, 1, 2, 3; indices k,l,m,n- ■ ■ — 1,2,3 and a,b,c,- ■ ■ — 0, 1, 2. The metric <? M „ has 

signature -\ and g is its determinant. Covariant derivatives are indicated by a D, partial derivatives by 

a d and covariant derivatives in a three subspace by V a or V&. The permutation symbol in 4 dimensions is 
e MI /p<r with eoi23 = 1 and in 3 dimensions e feim with £123 = 1. Finally rj^ pcT — ■ s f z ^ge ll , V p a . 



An extreme example of inertia due to gravitational waves alone is provided by Gowdy's [15] 
universe, a closed world that contains nothing except gravitational waves. One of our ultimate 
aims is to discuss the meaning of Mach's principle in General Relativity and particularly its 
application to such systems. However, such applications and the question as to what bearing 
Mach's principle has on the existence of dark energy are outside the scope of this paper. 

When we first discussed Mach's principle in FRW universes via perturbation theory [18] , 
we pointed out that our first order perturbation theory inevitably missed the second order 
contribution expected from gravitational waves (see also our recent comprehensive paper |3j ) . 
Here we show that there must be such a contribution. 

We consider waves in spaces that have a hypersurface orthogonal Killing vector along Oz. 
If such waves have axial symmetry their wave normals converge on the axis and they carry 
no angular momentum. To get truly rotating waves that carry angular momentum around 
the axis, we can not impose axial symmetry. As no exact non-axisymmetric rotating wave 
solutions are known, we work via perturbation theory. First we solve the flat space wave 
equation for linearized gravitational waves with symmetry along Oz, using cylindrical polar 
coordinates. Then we superpose these linear solutions to form a wave pulse that rotates 
around the axis. This pulse generalizes to no n- axial-symmetry the wave pulse found by 
Bonnor [7] and by Weber and Wheeler [30]. In the spacetime generated by such a pulse of 
linear waves, we then solve one of Einstein's equations to second order in the wave amplitude. 
Following [1 and [3], we see that the wave pulse generates an "effective source term" of second 
order in the wave amplitude. This means that the first order changes in the metric due to the 
wave can be neglected when calculating the extra perturbation due to these effective sources. 
(Keeping the changes could only lead to terms of third order in the wave amplitude.) Thus we 
have to calculate in flat space the effects of sources that are quadratic in the wave amplitude. 
Such perturbative techniques are well discussed in p3] ■ 

2 Calculations 

We consider as in [3] source free metrics with one hypersurface orthogonal Killing vector of 
spacelike translations £. Coordinates can be chosen so that z = x 3 and 

{x^} = {x a , x 3 } with a = 0,1, 2 that £ = {0,0, 0,1}, (2.1) 



and the metric can be written in the form introduced by Ashtekar, Bicak and Schmidt [T]: 

ds 2 = e- 2 ^g ab dx a dx b - e^dz 2 , (2.2) 

■0 and g ab are functions of the x c only. The source-free Einstein's equations R^ v = take 
then the following form in terms of the Ricci tensor 1Z ab of the 2+1-space da 2 = g ab dx a dx . 

R ab = => K ab = 2d a i>d b iP, (2.3) 

^33 = => g ab V a V b ^ = 0. (2.4) 



Because ip satisfies (2.4), the right-hand side of (2.3) automatically satisfies the Bianchi 
identities in 2 + 1 spacetime. To zero order g ab reduces to the flat spacetime metric and we 
choose cylindrical coordinates x° = £, x 1 = p, x 2 = (p, so, in the notations of [3] where a dot 
means a derivative with respect to t and a prime a derivative with respect to p, our linearized 
waves obey the wave equation 

$-il/'- V - \d 2 ^ = 0. (2.5) 

p p 2 * 

The solutions are expressed in Bessel functions J m {ujp): 

- l p = Ae ±imv±iujt J m (ujp). (2.6) 

We note that the special class of these waves with m = constitute the exact solutions 
mentioned in the Introduction. We follow [7] and [30] in choosing a superposition of such 
waves but in order to have rotating linearized waves we need m/0 and we choose waves of 
the form e " t ( rn{ P~^ t ) with mcu > so that the wave fronts rotate positively about the z axis. 
Our linear wave pulse is given by the superposition with a > being the effective duration 
of the pulse, 

roc 

ip = B (auj) m e- auJ e i( - mip - ujt h m (ujp)aduj + c.c, (2.7) 

./o 

B is a dimensionsless amplitude. It is convenient to define a = a{t) = a + it. Thanks to the 
Bateman Manuscript Project of Erdelyi et al |13j . formula 8.6.5, we find 

V; = Ba m+1 e im ^ • 2 m F(m + ^ • p m (a 2 + p 2 Y m ^ + c.c. (2.8) 

\/7T 



ij) may be written in terms of non-dimensional variables as in jj] , 

p=^ , t = -, and with 2 m T ^ m ^ 2 ' =(2m-l)!!, (2.9) 

a a V7r 



ip can also be written in real terms as follows: 



ip = 2B(2m - l)\\- 



p m cos [nop — (m + ^) x] 



ro \2 



[l + p 2 - t 2 ) z + At 2 



Km+i) 



where x = arctan 



It 



1 + p 2 - t 2 



. (2.10) 



So the phase of ip is imp — (m + 2) X- The p-dependence of this phase gives spirality to the 
wave and changes sign at t = 0. Figure 1 shows the profile of ip at different times for m = 10 
and ip = 0. 



To calculate the source term on the right-hand side of (2.3) we need derivatives of ex- 



pression (2.8). In particular the t,p component of the source will be given by 2ipd^ip. This 



product contains two terms that are independent of ip and two terms that vary like sin 2mp or 
cos 2mp. Those terms can not cause any rotation of the inertial frame on the axis. We there- 
fore concentrate on the terms that are independent of ip and denote them by the averaging 



symbol < >, implying an average over p. Defining j, we find from (|2.8|) and (|2.9|), 

j = -2 < ipd^ip > 
in which 



— m(2m + 1) [(2m - l)!!] 2 p 2m Q- m ^' 2 (^ + Ap] 
a \ou ) 



(2.11) 



ft and Q= (l+u-P) + At 



1+2 



(2.12) 

Before proceeding further we return to equations ( |2.3[ ), now regarding the right-hand 
side as a known second order source of perturbation of our now first order metric. As the 
source terms are already of second order, we may neglect the change of the metric from flat 
space in calculating their effects, as such changes could only give effects of third order. We 



can write the fully perturbed metric in the form (2.2) in which, see (5.4) in [3], 

g ab dx a dx h = e 27 (dt 2 - dp 2 ) - W 2 {dp - ujdt) 2 . 

The flat space equation for the axially symmetrical part of uj is then, cf. (5.10) in [3J, 

< ^02 >= — {ft <u»')' = 2 <ipd v ip>= -j. 
Ip 



(2.13) 



(2.14) 



We first integrate once with the boundary condition that <uj> is not singular on the axis: 



/T <uj> 



jdp\. 



Using the boundary condition that <u>^ at infinity we then find that 



<UJ>~- 



P2 \Jo 



jdpl dpi, 



(2.15) 



(2.16) 



where p\ and p\ are the dummy variables. Doing the outer integration by parts, we deduce 



<uj>-- 



1 
2m 



jdu + 



1 
2m 



jdu. 



(2.17) 



Substituting (2.11 ) into (2.17), we get, after integrating the Q m 3 ' 2 (dQ/du) terms by parts, 
2B 2 



<LO>~- 



in which 



-m[(2m- 1)!!] 



Im{u) = I tr 

'0 



/?? 



p 



I m _i + 2(2m + 1)— J m + (m - l)iJ m -i + 2(2m + l)riT„ 

M M 



l Q- m " 3 / 2 dM , H m (u) 



u 



-lQ-m-3/2 d 



U. 



(2.18) 



(2.19) 



The evaluation of <u» is thus reduced to the calculation of I m and H m for all m. This is 
done in the Appendix. 

3 Rotation of inertial frames at small and great distances 



On axis the first two terms in the large brackets of (2.18) are equal to zero. So using (A. 15) 
with b = 1 - i 2 and c = 1 + 1 2 we get from ( 2.18 ) 

B 2 (2m)! l + m(l + t 2 ) 



<UJ> -- 



(3.1) 



a 2 2 ™- 1 (1 + i 2 ) 2 ' 
This is greatest at t = exactly. Notice that there is no time lag between the wave arriving 
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Figure 1: This shows ij)/B, the wave profiles for m = 10 and tp = at 5 different times at adapted 
scales. On the left-hand side ip/B is shown as it circulates inwards, from right to left at times 
i = —10,-6, and —3. The right-hand side shows it as it circulates at closest approach to the axis 
t = 0. Notice there is very little closer than p ~ 0.4 even at t — 0. Superposed on the right-hand 
side is the profile as the wave starts to retreat at i = 3. Thereafter the profiles arc the same as on 
the left-hand side of the figure but in reverse order. At all times the wave profile centers around 
p= Vl + i 2 - 



closest to the axis - most of the energy never gets nearer than p ~ 0.4 (see Figure 1 at 



t = 0) - and the resulting rotation of inertial frames. Just as in the work of Lindblom 
and Brill [16J on inertia induction from falling rotating spheres, the non-local effect appears 
instantaneously when suitable Machian coordinates (gauge) are chosen. This is a consequence 
of Mach's principle being embodied in the "instantaneous" constraint equations of General 
Relativity. In electricity, the presence of a charge within a sphere at time t can be ascertained 
"instantaneously" by measuring the flux of its field out of a surrounding surface. You do not 
have to wait for the field to come out from the charge - it is there already. 

Returning to our inertial frame on the axis, a rod lying peacefully at the origin, initially 
oriented at ip = cpo, will remain in the flat space near there and feel no torque, but, due to 
the rotation of the inertial frame it will point towards tp = <fi(t) where 

f 



<t>{t) =0(0) + 



* , , , ,(2m + l)! 

<uj> dt = ( f ) (0)+B 2 y 92m ; 



TT 

arctant -\ \- - — 

2 (2m + l)(l+t 2 ) 



When the whole pulse has passed it ends up oriented towards 



(oo) = 0(0) + B 



2 (2m + l) 

2 2m 



TT. 



For symmetry we choose 



so that = when t = and 



0(0) 



-B 



(2m + 1)! vr 
2~2™ ' 2 ' 



6{t) = B 



2 (2 771+1]! 

22m 



arctan t + 



(2m + l)(l + t 2 ) 



(3.2) 



(3.3) 



(3.4) 



(3.5) 



Figure 2 shows 0(t) for t > 0. It is antisymmeric about t = 0. 



Far from the axis we have from (2.18) and (A. 9) 

B 2 m(ml)(2m - 1)11 



1 



/5> 1. 



(3.6) 



a 2™- 1 p 2 

As expected this is constant in time for p large, but falls off like pT 2 as befits a cylindri- 
cal system with a total angular momentum per unit height, cf. [1], 

tx 2 m(m\)(2m — 1)!! 



J 



-aB< 



)m—2 



(3.7) 



For general t and p, <u>> is given in (2.18) with I m and H m from (A. 6) and (A. 12). But 
for p 2 <C 1 + i 2 we need to show that our space is almost flat so that < u > can be properly 
interpreted as a rotation of the inertial frame. We therefore evaluate <cu> not just at p = 0, 




Figure 2: The angle <fi(i)/B 2 through which the inertial frame at the origin turns as a function of 
time for m = 10. The initial angle before the wave comes in, is chosen as — 0o so that cf> is zero at 
t = 0. This produces anti-symmetry about t = 0. w(i,0) is given by the gradient of the curve. It is 
greatest at t = 0. 



but in this larger region, and show that at each time < w > varies weakly with p there. Indeed 
from ( |Al7) and flAT8| ) 



<a;>~<w> 



(2m -1)!! 



(2m - 5) + (2m + 7)£ 2 



. (3.8) 



(m + 1)! (l + t2)m+i (m + 1) +mt 2 V 1 + * 2 . 

Figures 3 and 4 show j, k times the angular momentum density, and the resulting <u(t, p) >. 



The fact that the source terms in (2.3) all vary as p shows that the other terms in the 



metric will vary similarly, so wherever [p 2 /(l +P)] m is negligible our second order metric will 
take the form 



ds 2 = e 2 ^ 



dt 2 - dp - p (dip - < uj > dty 



e^dz 2 



p 2 <.{l + t 2 ) (3.9) 



with t/j given by (2.10). If we now go to rotating axes by writing (p = ip — 4>(t), then our 
metric becomes 



ds 2 = e 2 ^ (dt 2 - dp 2 - p 2 d(p 2 ) - e^dz 2 



p 2 < (1 + if) 



(3.10) 



which is flat in the reduced space. Furthermore the wave itself never penetrates significantly 



within p/V 1 + t 2 ~ 0.4. 
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Figure 3: The profile of the angular momentum density at t ~ 0, 1 and 2. Again the drawing is for 
m= 10. 

4 Angular momentum transport by gravitational torques 

Those who have studied angular momentum transport in spiral galaxies, know that the 
classical gravitational stress tensor is [T7] , [25] 



a kt = -(2d k tpdiijj - S k i y~) \d m i>f 



(4.1) 



where here ip is the classical gravitational potential. The gravitational couple transfering 
angular momentum outwards through a cylinder of radius p is 



C„ 



£3kiX (a m dS m ), 



(4.2) 



where dS is the outward pointing surface element and X the radius vector from the axis. 
Evidently this reduces to 



c, 



dptpd^tppdipdz. 



(4.3) 



For these gravitational stresses to carry angular momentum outwards there must be a positive 
correlation between d p ip and d^ip averaged over such a cylinder. Such a correlation gives a 
trailing sense of spirality to contours of constant ip in the sense that the outer parts of a spiral 
galaxy trail that of the inner parts in the sense defined by the rotation. It is no different 
here; for angular momentum transport outward, the spiral formed by contours of ip must 
trail, but contrarywise when they form a leading spiral with the outside further advanced 
than the inside, the angular momentum is transported inward. Thus our knowledge of spiral 
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Figure 4: The profile <w> as a function of p at selected times i = 0, 0.5, 1, and 2. Notice the still 
flat uniformly rotating central region when t = 0. 

structure in galaxies enabled us to predict that there would be leading spirals when the wave 
was going in and therefore transporting its angular momentum inward, but trailing spirals 
when our gravitational wave was receding from the axis. At t = there should be no angular 
momentum transport so the contours of ip should form a rotating cartwheel structure with 
no spirality. These expectations were beautifully confirmed when we plotted the contours of 
ip using Mathematica. The waves always rotate counterclockwise in Figure 5. They start 
as quite tightly wrapped leading spirals; as they move inwards they keep rotating but the 
wrapping becomes less tight. At t = they unwrap to form a rotating cartwheel structure 
and thereafter they move outward as trailing spirals carrying angular momentum with them. 
It is when they are closest in that their rotation causes the greatest rotation of the inertial 
frame in the central region, though the waves themselves never reach much further in than 
p = 1. Nevertheless the second order effect of the angular momentum causes the rotation of 
the inertial frame within. This is illustrated by the direction of the central arrow that shows 
the orientation of a rod initially at rest at the origin. Because this effect is second order, 
the rotation seen in the figure is dependent on the amplitude of the wave, whereas the spiral 
shapes are not. 
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5 Energy and angular momentum 

In Section 3 we used the asymptotic form of < uj > to deduce the total angular momentum per 
unit height but we could also calculate it via the angular momentum density j/k 

1 f°° 4vr „, ., ,,,., r T , nT9t , vr nm(m!)(2m-l)!! 



1 f 00 47T 7T 

J = - j2irpdp = — aB 2 m[(2m-l)ll} 2 \ml m -i + (2m + l)Pj m l = -aB 2 

Kj K K 



2m— 3 

(5.1) 



Any quantity ip that obeys Oip = automatically conserves both J and the quantity 

1 



The identity 



(^ z + \Vip\ 2 )dV. (5.2) 

K ' 



-d t {i> 2 + |VVf) = -V • (2^V^) (5.3) 

K K 



follows directly from the wave equation and can be interpreted as a conservation law with 
-ipVip ■ dS being the flux through an element of the closed surface surounding a volume V. 
We notice that the flux vector now directed through a surface ip =const will on multiplication 
by p give the angular momentum density j/k. Integrating -(ip 2 + IV^I 2 ) over a volume of 
unit height we find the grand total and using the formulas of the Appendix, 



OO fllX 



1 f.l.2 , l^/.l2^ J ,„„ J „_^^2^K2m + l)!! 

^0 



{ip 2 + \V^\ 2 )dppdp = -B 2 ■' r \ (5.4) 

K K 2 m l 



which is independent of time as expected. We may identify this conserved quantity with 
energy per unit length via the following argument. The conservation laws Vfe^ = 0> see K> 
in the reduced 2 + 1 space follow from its contracted Bianchi identities of which the x° = t 
one reads 

l - [(pT °y + (pTo 1 )' + <y pT 2 )] = 0. (5.5) 

Ordinary derivatives in flat space have replaced covariant ones because T£ is already of second 
order in the wave amplitude: 

Kl ah = TZ ab - \g ah ll = 2d a ipd h ip - g ab (V 2 - |VVf ) , (5.6) 

where Vip is the normal flat space gradiant in p and tp, 

|VVf = V' 2 + 4(<V/f- (5-7) 



From (5.6) we see that 

kT ° = i; 2 + |VVf (5.8) 

11 



which is the density of our conserved quantity. The total energy per unit length as measured 
from the conicity at infinity, 700, was shown to be Ik^^Jk by pQ. 



In [3], eq. (5.31), we show that Its^^Ik is given by the left-hand side of (5.4). This 
completes the identification of our conserved quantity with the energy per unit height. It 
also shows that the wave produces a long range effect on 7 far from the region p ~ t where its 
main amplitude lies. It is one of the strange properties of cylindrical General Relativity that 
this conicity is independent of p outside the wave. This has the effect that while particle 
trajectories are globally bent by the defect angle nevertheless there are no orbits along finite 
sections of circles (about the axis) caused by the long distance gravity of the cylindrical wave. 
Rather the geodesies are like the geodesies on a cone all of which go to infinity. It would 
be interesting to calculate the gravity due to a more realistic three dimensional wave pulse 
for which the odd behaviour at infinity found in 2+1 dimensions would not obscure the 
interpretation pHT| . 
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A Appendix: Evaluation of the Integrals I m and H m 

ru poo 

I m = u m Q- m -^ 2 du ; H m = / u^Q-^^du. (A.l) 

JO Ju 

Q = u 2 + 2bu + c 2 , where b = l-P and c = l + t 2 , (A.2) 
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Now differentiate m times with respect to 6: 



O- 2 ,7„ ) = ^Q- 1/2 - ^ = *(&,«)■ (A.4) 
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So 
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(-l) m d m X(6,u) 



(2m + 1)!! 56 
where in this expression c = 1 + t 2 and, after differentiation b = 1 — i 2 . Now from (A.4) 
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Differentiating with respect to c 2 and multiplying by —2/3: 



/ Q~ 5 / 2 di 



1 



3(c 2 - 6 2 ) 2 



2 f,2\ 



{u + b){2Q + c 2 - b 

Q3/2 



y(M). 



Differentiating m — 1 times with respect to 6 we find 
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Hence, 
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,ra- 1 



(A.12) 



again b must be set equal to 1 — t after the differentiation. From (A. 10) we see that 
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Whereas (A. 6) and (A. 12) give / m and H m for all values of u, (A. 8) and (A.15) give their 



values at and oo, nevertheless it is useful to view their behaviours at small and at large 
u explicitly. In particular for p 2 <^ ^(1 -\- i 2 ), the region encompassed by the wave but not 
'feeling' it, we may expand Q~ m ~ 3 > 2 in the form 
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Finally when b = c, Q becomes a perfect square so the integrals are easier; since both b and 
c are then 1 we express the answer that way: 
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Figure 5: This shows the m = 10 wave which always rotates anticlockwise. As it comes inwards 
(Figure 5a) at t = —1 it is in the form of a leading spiral with the outer parts of the arms ahead of 
the central parts. By t = — 1 (Figure 5b) the spiral has started to open. By t = (fig. 5c) the central 
parts have caught up and the spiral has changed to a cartwheel structure but rotation keeps it beyond 
p ~ 0.4. By t = 1 (Figure 5d) the spiral has become trailing as befits a wave that now feeds angular 
momentum outwards. By i = 2 (Figure 5e) the spiral becomes tighter and the flat central cylinder 
becomes larger. We show i — 10 (Figure 5f ) at a small scale but notice the beautiful tight wrapping 
of the narrow arms. Also notice the opposite spirality of the conjugate pairs t = ±2 and t = ±1. The 
figures encompass a radius p ~ 7 (p ~ 17 for i = 10). The height of ip/B reduced by a factor of 10 -4 



